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Abstract 

Optimal truncations of asymptotic expansions are known to yield approximations 
to adiabatic quantum evolutions that are accurate up to exponentially small errors. In 
this paper, we rigorously determine the leading order non-adiabatic corrections to these 
approximations for a particular family of two-level analytic Hamiltonian functions. Our 
results capture the time development of the exponentially small transition that takes 
place between optimal states by means of a particular switching function. Our results 
confirm the physics predictions of Sir Michael Berry in the sense that the switching 
function for this family of Hamiltonians has the form that he argues is universal. 
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1 Introduction 



The adiabatic approximation in quantum mechanics asymptotically describes solutions to 
the time dependent Schrodinger equation when the Hamiltonian of the system is a slowly 
varying function time. After a rescaling of the time variable, the adiabatic approximation 
describes the small e behavior of solutions to the Schrodinger equation 

ie^ = H(t)^. (1.1) 

In the simplest non-trivial situation, {H(t)} t£ ^i is a family of 2 x 2 Hermitian matrices 
that depends smoothly on t, and whose eigenvalues Exit) and E 2 (t) are separated by a 
minimal gap E 2 {t) — Ei(t) > g > for all t e JR. 

To discuss scattering transition amplitudes, we also assume that H(t) approaches limits 
as t tends to plus or minus infinity. We let $?(£), j — 1, 2 be smooth normalized instantaneous 
eigenstates associated with Ej(t), respectively. Then the transition amplitude A(e) across 
the gap between the asymptotic eigenstates is defined as 

.4(e) = lim |($ 2 (ti), U e {ti, *o)$i(to))|, (1-2) 

[Q — * — oo 
t\ — *- + oo 

where U t (ti, to) denotes the evolution operator corresponding to (jl.ljl . The adiabatic the- 
orem of quantum mechanics, 6\, asserts that A(e) = 0(e), so the transition probability 
A(e) 2 is of order 0(e 2 ). 

If the Hamiltonian is an analytic function of time, the transition amplitude is much 
smaller. Long ago, Zener p>3; considered a specific real symmetric two-level system, that 
had an exponentially small transition A(e) as e — > 0. Generalizations of Zener's result 
to analytic real symmetric two-level Hamiltonians with a gaps in their spectra were then 
proposed in the physics literature. Formulas for A(e) of the form 

A(e) ~ e~ 7/e , as e -> 0, (1.3) 

with 7 > that applied more generally were obtained, e.g., in [21], 0, jUj. The decay rate 
7 was essentially determined by complex crossing points, i.e., the points in the complex t- 
plane where the analytic continuations of the eigenvalues coincided. Later on, papers [3] and 
[TS] recognized independently that non-trivial prefactors G, with \G\ ^ 1 could be present 
for general Hermitian two-level Hamiltonians to yield the general formula 

A{e) ~ |G| e" 7/£ , as e -> 0. (1.4) 

Also, 5J and [T3] pointed out independently that certain complex degeneracies could lead to 
the same formula in the real symmetric case. Formulas (jl.3j) and (jl.4j) are variants of the 
well-known Landau-Zener formula that has been widely used in many areas of atomic and 
molecular physics. 

The goal of this paper is to obtain more precise results for a family of two-level systems. 
For all times t, we construct an approximate solution to (jl.l|) that is accurate up to errors 
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of order exp(— j/e) e M , for some (i > 0. This captures the transition process which is of order 
exp(-7/e). 



Explicitly, let E > and 5 > 0, and consider the Hamiltonian function 

Hit) = E ( J * ) (1.5) 

whose eigenvalues are ±E/2 for every t. This Hamiltonian can be viewed as the familiar 
Landau-Zener Hamiltonian modified to keep its eigenvalues constant. In our case, the notion 
of eigenvalue crossing point is replaced by the singularities of the Hamiltonian itself at 
t = H 5. These points govern the transitions between the two levels. Let and $ 2 (^) 

be smooth, normalized real eigenvectors corresponding to —E/2 and E/2 respectively. 

Theorem 1.1 Let H(t) be given by M.b}) and let < /i < 1/2. Then: 

1) There exist vectors Xi{ e -> an d X2{t, t) that satisfy the Schrddinger equation M-l)) up 
to errors of order e~ ES ^ e and correspond to the eigenstates and $ 2 (t) in the sense that 

lim | t))\ = 1+ 0(e-^/^). (1.6) 

Moreover, the set {Xj( e > t)}j=i,2 is orthonormal up to errors of order e~ ES ^ e e M . 

2) The Schrddinger equation has solutions \Pj(e, t), j = 1, 2 such that uniformly int G M 
as e — > 0, 

^( e ,t) = Xl (e,t) - v^e^^i jerffy^tj + l| X2 (e, f) + 0(e~^M 



and 



* 2 (e,f) = X2 (e,t) + vSe"*'* i j erf ^ tj + lj Xl (e, t) + 0(e- £ ^) 



Remarks 

0. Recall that the function erf is defined by 

erf(x) = — f e~ y2 dy = — [* e~ y2 dy - 1 e [-1, 1]. (1.7) 

^/7T JO y7T J-oo 

1. The vectors Xj( e ! 0> j 1 = 1>2, are constructed as approximate solutions to (jl.lj) obtained 
by means of optimal truncation of asymptotic expansions of actual solutions. As t — > — oo 
they are asymptotic to the instantaneous eigenvectors $j(t) of H(t), up to a phase. We call 
them optimal adiabatic states. See (|4.1Uj) . (|4.11|) and (|6.9|) . 

2. The transition mechanism between optimal adiabatic states takes place from the value 
zero to the value \/2 e~ ES ^ e in a smooth monotonic way described by the switching function 
(erf+l)/2, on a time scale of order ^fe. By contrast, the transition between instantaneous 
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eigenstates of the Hamiltonian displays oscillations of order e for any finite time to eventually 
reach its exponentially small value only at t — oo. In that sense also, the vectors Xj{ e i t) are 
optimal. 

3. As the optimal adiabatic states and eigenstates essentially coincide at t = ±00, the 
transition amplitude equals A(e) ~ \/2e~ E& ^, up to errors of order e~ ES ' e e (Jl . 

4. The parameters E and 5 play somewhat different roles. If we fix E and decrease S, the 
singularity of Hit) approaches the real axis. The transition amplitude \/2e~ ES ^ e increases 
whereas the time it takes to accomplish the transition decreases as \jlbejE. If instead we fix 
5 and let E decrease, the gap between the eigenvalues decreases. The transition amplitude 
increases as well, whereas the typical time of the transition now increases. 

5. Our results allow us to control the evolution operator U e (t, s) associated with up to 
errors of order e~ E5 ^ € e M , for any time interval [s, t]. 

6. Further comments concerning the relevance of the Hamiltonian (jl.5|) are presented at the 
end of this section. 

We now put our results in perspective by describing previous work on exponential asymp- 
totics for the adiabatic approximation. 

Rigorous computations of the behavior of the exponentially small quantity -4(e) for two- 
level systems, or generalizations of this typical setting, were provided relatively late in |18j . 
[2U] . [To] . [To] . [2Tj . [T7j . [2B] . Although we shall not use the technique in this paper, let us 
briefly describe the mechanism that is typically used to get the asymptotics leading to the 
exponentially small quantity -4(e) in these papers. It involves deforming integration paths 
from the real axis to the complex plane t-plane until they reach a (non-real) crossing point. 
Crossing points provide singularities where significant transitions take place, but their lying 
away in the complex plane makes these transitions exponentially small due to the presence of 
dynamical phases exp(— i / * Ej(s) ds/e) whose exponents acquire a non-zero real part along 
the path. With this approach, the link with the initial problem posed on the real axis is 
possible only at infinity. One does not learn about the dynamics of the transition. We note 
that in a more general framework where -4(e) represents the transition between two isolated 
bands of the spectrum for general (unbounded) analytic Hamiltonians, exponential bounds 
on -4(e) were obtained by suitable adaptations of this method in [15] . [T4] . See also [HH] 
for similar results using the pseudo-differential operator machinery. 

The other successful method used to construct precise approximations of solutions to (jl.lj) 
uses optimal truncation of asymptotic expansions. With sufficiently sharp estimates of the 
errors, one can prove exponential accuracy. Under appropriate analyticity assumptions, one 
typically proves that the error committed by retaining n terms in the asymptotic expansion 
is of order n\ e n . This error is minimized by choosing n ~ 1/e. By virtue of Stirling's formula, 
it is of order exp(—l/e)/y/e as e — > 0. The optimal truncation method was first used for 
the adiabatic approximation by Berry in [3|. He constructed approximate solutions to (jl.lj) 
and gave heuristic arguments concerning their exponential accuracy and the determination of 
the exponentially small transition mechanism between asymptotic eigenstates. In particular, 
the switching function (erf+l)/2 in Theorem 1.1 first appeared in Berry further claimed 
that this function was universal, i.e., the time development of non-adiabatic transitions in all 
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systems governed by one crossing point (and its conjugate) were described by this switching 
function. His formal arguments were supported by beautiful numerical investigations of Lim 
and Berry [25] . 

Berry's paper [Hj is the main inspiration for the present work. 

Mathematically rigorous exponential bounds for solutions to (jl.ljl (and of -4(e)) using 
optimal truncation were first obtained in the general situations by Nenciu in |3T] . They were 
refined later in [20]. See ^T] for an elementary derivation of such results. 

Although these rigorous results prove the exponential accuracy of the optimal truncation 
technique, their estimates are not accurate enough to capture the exponentially smaller non- 
adiabatic transitions. An exponentially small bound on -4(e) is an easy corollary, but the 
estimates do not provide the asymptotic leading term (jl.4|) for -4(e). 

We refer the reader to two fairly recent reviews, [22] and PP, for more details and many 
other aspects of the adiabatic approximation in quantum mechanics. 

We also note that in the broader context of singular perturbations of linear ODE's, 
Theorem 1.1 can be interpreted as the smooth crossing of a Stokes line that emanates from 
some eigenvalue crossing point or singularity. See e.g., [2], |2ZI, [III [Hj and references 
therein. However, from our perspective, in all rigorous work that has dealt with such issues, 
the crossing of a Stokes line is performed on a very small circle around the point responsible 
for the transition. In this paper all estimates are performed on the real axis, and the two 
conjugate points responsible for the transition are fixed and away from the real axis. 

Another angle of attack for such problems uses Borel summation ideas. This is done for 
certain singularly perturbed ODE's e.g., in The method consists of writing the solution 

as a Laplace transform evaluated at 1/e, i.e., as / F(t, p) e P//e dp. One then derives a PDE 



for F. The PDE is roughly what one obtains from the original equation by replacing 1/e by 
the symbol d/dp. To get exponentially precise information about the solution to the original 
problem, it is enough to study the location and nature of complex singularities of F. 

We tried to implement the Borel summation technique for our problem, but failed to 
obtain sufficiently detailed information on the nature of the singularities. 

To the best of our knowledge, no rigorous results that address the issue we describe in 
Theorem 1.1 are available in the literature. 

Before we turn to the proof of Theorem 1.1, let us briefly discuss our choice of Hamiltonian 
(II. 5|) . This choice belongs to the family of real-symmetric time-dependent 2x2 Hamiltonians 
with non-degenerate eigenvalues E 2 (t) > E\(t). For any member of that family, we can 
assume without loss of generality that E 1 (t) = —E 2 (t), because we can subtract a time- 
dependent multiple of the identity from H(t), which only changes the solutions by a trivial 
phase. If we change the time variable from t to t' — 2 Jq E 2 (s) ds/E and drop the prime on 
t' , we obtain a Schrodinger equation (jl.ljl with a new Hamiltonian h(t) of the form 




(1.8) 
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whose eigenvalues are ±E/2 for every t. The angle a is given by some function of time. We 
choose 

^K*)/2)Y , . ( oos(a(0/2) A 



= I cos(a((t)/2) J ' and $2(t) " WnH(0/2) ; 

as the normalized real eigenvectors of Then the coupling f(t) that drives transitions 
between the instantaneous eigenvectors is given by 

f{t) = (* 2 (f),*i(t)> = -<*i(*), = 
Our choice of Hamiltonian (|1.5|) corresponds to 

/(*) = 2 ^2 + p\ ° a W = ~ ^ arctan(t/<5), 

where 5 > is a parameter monitoring the strength of the coupling. This choice of coupling 
presents the simplest non-trivial singularities in the complex t-plane. 

From another point of view, our choice of Hamiltonian is motivated by the Landau- 

Zener Hamilton, ( \ ) , which has the local strnctnre of a generic avoided crossing 

[TU] . Physically, one expects the transition to take place in a neighborhood of t = 0, and 
one expects that only the form of the Hamiltonian for small t should determine the transi- 
tion dynamics. To first order near t — 0, our Hamiltonian agrees with the Landau-Zener 
Hamiltonian. 

However, quantum transitions have a more global character. The presence of the square 
root factor in (jl.5j) . gives rise to the nontrivial prefactor v^2 in the transition amplitude. 
Indeed, it is shown in ^H] that a change of variable allows one to transform equation (jl.lj) 
with Hamiltonian (|1.5J) into an equivalent Schrodinger equation driven by a Hamiltonian 
that behaves locally near t = as a Landau-Zener type Hamiltonian, but with a non- 
generic complex crossing point. This non-generic structure is responsible for the nontrivial 
prefactor \/2 in the transition amplitude. It follows from fH], Section 4, that the leading 
order transition amplitude for our Hamiltonian is 

A(e) ~ V2e~ ES/e . (1.9) 



The rest of the paper is organized as follows. In the next section we develop perturbation 
expansions for solutions to the time-dependent Schrodinger equation (jl.lj) with Hamiltonian 
(jl.5j) . In Section |21 we analyze the behavior of the high order terms of this expansion that are 
required for precise optimal truncation. In Section 0] we study the error term obtained from 
optimal truncation and define optimal adiabatic states. Section IS] is devoted to the study of 
two integrals that arise in the error term and give rise to the switching function. Theorem 
1.1 is then proven in Section |UJ 

Acknowledgements George Hagedorn wishes to thank the Institut Fourier of the Univer- 
site de Grenoble I for its kind hospitality and support. Alain Joye wishes to thank Virgina 
Tech for its kind hospitality and the NSF for travel support. We also wish to thank Ovidiu 
Costin for many useful discussions about this problem. 
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2 The Formal Perturbation Expansion 



We start by converting our time-dependent Schrodinger equation (jl.lj) with Hamiltonian 
()1.5|) into a parameter free equation. 

Changing variables from t to s with s — t/8, we obtain 

. (ty E5 < 1 s 

ie ds ~ 2VTT^ V s _1 

Thus, without loss of generality, by changing e into e' — e/ we can study the parameter 
free model with 5 = = 1 in (jl.lj) . We drop the prime, keep e in the notation, and consider 
the coupling 

/W = J^Wy (2.D 




We now develop a formal asymptotic expansion to solutions of fll.ljl . 

We concentrate on constructing a formal perturbation expansion of the solution to (jl.lj) 
that corresponds to the negative eigenvalue —1/2 for small e. We make the unusual ansatz 
that (jl.lj) has a formal solution of the form 

^(e, t) = e it/{2e) efo f( s )^,s)ds ( $ ^ + g ^ t) ^ j ^ ^ 

oo 

where #(e, t) = ^ 9 3 We^ 
i=i 

Remarks 

1. We arrived at this ansatz by attempting a formal solution of the form 

,j,(e,t) = e^e 2 ^ ($ x (f) + p(e,t) $ 2 (t) ) , 

with z(e, t) — Zj(t)ei. We then realized that this required z(e, t) = f(s)g(e, s)ds. 

3=1 

2. There are more standard ansatze for the perturbation expansion [3j ITU ITo! ITS! I19j !2Uj 

122 12S 123 123 EOl ED E3, but we were unable to get sufficient control of their n th terms to 
prove the estimates that we required. 

3. If we do not expand g(e, t), then we find that it must satisfy 

ieg'(e, t) = g(e, t) - ief(t) (l + g(e, t) 2 ) . 
However, we will not use this equation. 

4. For normalization purposes, we later consider (j2.2j) with JL^ f(s)g(e, s) ds in the expo- 
nent, instead of Jq f(s)g(e, s) ds. 

5. When seeking a solution that corresponds to the positive eigenvalue for small e, one 
makes the similar ansatz 

^( e ,t) = e - it/{2e) e ~ & f{s) ~ 9(t ' s)ds ($ 2 (t) + £(e, t)*i(t)), (2.3) 
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where g(e, t) — gj(t) e J satisfies 
i=i 

ze#'(e, t) = -g(e, t) + ief(t) (l + g(e, t) 2 ) . 
Hence, for any j G N, we have gj(t) = gj(—t). 

We substitute ()2.2|) into (jl.lj) and formally solve the resulting equation order by order in 
powers of e. 

First Order: The terms of order e require 

9i(t) = if(t). (2.4) 
Second Order: The terms of order e 2 require 

g 2 (t) = ig[{t). (2.5) 
Third and Higher Order: The terms of order e n+1 for n > 2 require 

9n+i(t) = i (g' n (t) + f(t) E fcW^Wj. ( 2 -6) 

Using ()2.1|) for the coupling /, and an easy induction using partial fractions decomposi- 
tions, we see that g n {t) can be written as 

2n 

9n(t) = E c ">; e i(*)> where e 2j _!(t) = (l + zt)- J ', e 2i (t) = (1 (2.7) 
i=i 

Thus, for each n, we can associate g n with a unique element of Z 1 , the space of absolutely 
summable sequences. 

Following the intuition of Michael Berry jHOH], we isolate the highest order poles of g n {t) 
by decomposing g n (t) = G n (t) + h n (t), where G n (t) = c nt2n ~i e 2n -i(t) + c nM e 2 n{t)- 

Proof of our results now depends on an analysis of the behavior of G n {t) and h n (t) for 
large n. 

3 Analysis of the Recurrence Relation 

The main goals of this section are summarized in the following proposition. 

Proposition 3.1 There exists C, such that for all real t and each n > 1, 

\G n (t)\ < (n-l)\, (3.1) 

\G' n (t)\ < n\, (3.2) 

\h n {t)\ < C(n-2)! log(n-2), (3.3) 

K(t)\ < C(n-l)!log(n-2). (3.4) 
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We prove this result using the recurrence formulas ()2.4|) - (j2.6j) and the following amazing 
fact that we use to control the nonlinear terms in the recurrence relation. 

fc+m+l 

Lemma 3.1 For each k andm, e k {t) e m {t) = ^ 4,mj e j{t)- 

i=i 

k+m+l 

Every d k ^ m j is a non-negative real number, and V] 4,mj = 1- 

j=i 

Remark This lemma implies that we have a Banach algebra structure on I 1 , where the 
product of {a n } and {b n } is determined by formally multiplying J2 n a n G n (t) times J2 n b n e n (t) 
and then taking the coordinates of the result in the {ej(t)} basis. 

Proof of Lemma 13.11 We first remark that by keeping track of the orders of the poles, it 
is easy to see that j < k + m + 1 in the sums in the lemma. 

Next, we observe that there are many trivial situations. If k and m are both odd, then 
ejfc(t) e m (t) = e k+m+1 (t). If k and m are both even, then e k (t) e m (t) = e k+m (t). 

Thus, the only non-trivial cases are when one is odd and the other is even. To prove 
these cases, we do inductions on odd k and even m. 

For k = 1 and m = 2, the lemma follows immediately from 

(1 + it)" 1 (1 - it)" 1 = \{l + tt)- 1 + i (1 - it)" 1 . (3.5) 

Next, we fix m = 2 and assume inductively that the lemma has been proven for k = 
2K — 1. Then using (|3.5jl again, we have 

e k+2 (t)e m (t) = (I + it)-*- 1 (I -it)" 1 



;i + it) 



-K 



= ^ [ e *+2(*) +e k (t)e 2 (t)} (3.6) 

The result now follows from our induction hypothesis. 
Thus, the lemma is true for all k and m = 2. 

Finally, we fix an odd k, and assume inductively that the lemma has been proven for this 
k and an even m. Then 

e k (t) e m+2 (t) = {e k (t)e m (t)) e 2 (t) 

= d k ,ra,je j (t)e 2 (t) 
j 

= 22 dk, m jzl dj,2j'ej'(t). 
j j' 

The lemma now follows since X^' dj t2 ,j' = 1, J2j d kjm j = 1, and all the d's that occur 
here are non-negative. | 

We also need the following technical result. 
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Lemma 3.2 For each n > 1, we have 

{n-jy.f 



^v^MUL < |, (3.7) 



n-1 



j=0 

n-1 



n\ 




(n - j) 


! j\ 


n\ 




(n - j) 


! j\ 


n\ 



- 3' anc? 



2 
3' 



Proof The first inequality trivially implies the other two. 

We observe by direct computation that the result is true for the first few values of n, and 
that the sum equals 8/3 when n = 3 and n = 4. 

For n > 5, we separate the first two terms and last two terms to see that the sum equals 

2 + - + 



n £2 ra! 

The largest terms in the sum over j in this expression come from j = 2 and j = n — 2. 

2 

Those terms equal — -, and there are (n — 3) terms. Thus, the left hand side of (|3.7|) 

n(n — 1) 

is bounded by 

2 2(ra -3) 2 1 8 

2 + - + -) { <2 + - + -<-. 

n n(n — 1) 5 5 3 

This last step relies on the observation that ^(n-i) takes the value 1/5 when n = 5 and n — 6, 
and that it is decreasing for n > 6. | 

For any = J2j Uj e jif)^ with G I 1 , we define ||y|| = J2j \Uj\- We note that for 
t e M, \y(t)\ < \\y\\. Since G n is obtained from g n by dropping components in the Bj(t) 
basis, we note that \\G n \\ < \\g n \\- Thus, the following lemma implies (|3.1|) and (|3.2|) since 



[l±it)- j = Tij (l±it)- j ~ 1 . 



Lemma 3.3 \\g n \\ < (n— 1)!. 



Proof We prove that the sequence a n = \\g n \\/(n — 1)! is bounded above by 1. 
By (J2.6|) . Lemmas 13. II and 13.21 we see that n > 2 implies 

O-n+1 < a n + 7T— : rr- (3i 

3n(n — 1) 



From (|2.4|) and (|2.5j) we have ai = a 2 = 1/2. By explicit computation, we observe that 

^3 = I2 - 1 _ 9' alld ° 4 = 384 - 1 ~ 12- 
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The lemma now follows by induction (starting at n = 4) and the following statement: 
If n > 4, a n _i < 1 - and a n < 1 - ^, then a n+ i < 1 ; 

To prove this statement, we use (j3.8|) to see that for n > 4 



3(n+l)' 



4 4 

OVt+1 < 1 - h 



1 - 



= 1 
< 1 



3n 3n(n — 1) ^ 3(n — 1} 
4 8(3n 2 + lOn - 29) 



3(n + l) 27(n+ l)n(n - l) 3 
4 _ 



3(n + l)' 

Lemma 13.11 is now proven by the comments before Lemma 13.31 and the following lemma. 

Lemma 3.4 There exists C , such that \\h n \\ < C (n — 2)! log(n — 2) and 
\\K\\ < C(n-l)! log(n-2). 



Proof The first estimate implies the second since h n is in the span of ej(t) with j < 2n — 2. 

Define b n = \\h n \\/((n — 2)!). Since h n is obtained from g n by dropping components in 
the ej(t) basis, we have \\h n \\ < \\g n \\- Thus, by Lemma 1331 b n < n — 1. 

We rewrite (|2.6|) . using g n = G n + h n : 

dG dh n_1 n_1 

C n +i + /i n +i — i \ ^ + + / G n _j + 2/ + f hj h n _j 



i=i 



i=i 



We then drop the e 2 n+i(t) and e 2n +2(^) components of this expression to obtain an expression 

dG 

for h n+ i. This involves dropping the entire term i — - — , as well as parts of other terms. Since 

dt 

the norm decreases when we drop components, we see that 



\h 



71+1 | 





dh n 








+ 




dt 



n-l 
3=1 



+ 2 



n-l 

i=i 



n-l 



dt 



Since 7t n is in the span of {ej(t)} for j < 2n — 2, 
\\Gn\\ < (n — l)\, h± = h,2 = 0, Lemmas 13. lj and 13 .2\ we have 

5/3 1/3 



< {n-l)\\h n \\. Thus, by ll/H < 1/2, 



4 /3 

n +l < n H 7 + 



&n-l + 



n-l (n-l)(n-2) (n - l)(n - 2)(n - 3) 

Since we already have b n < n — 1 and b\ = hi = 0, this implies 

4/3 5/3 1/3 , 10/3 
b n +i < K + — + — '— + — ^— = 6 n + - — '—, 



n — 1 n — 1 n — 1 



(n-l) 
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for n > 2. Thus, 



n— 1 n— 1 



6n = E (fcfi = Er4 < T (1 + log(n-2)). 

This implies the lemma and completes the proof of Proposition 13.11 | 

The functions ej(t) are in L l (M) for j > 3. The functions exit) and e2(i) are not in 
L x (iR), but ei(t) + e2(t) is. The following lemma facilitates getting L 1 information about 
the functions g n . 

Lemma 3.5 For every n, the e\(t) and ez(t) coefficients in g n (t) are equal. Furthermore, 
the L^-R) norm of g n is bounded by n \\g n \\. These results are also true for h n . 



Proof We prove the first statement by induction on n. It is clearly true for n = 1. Suppose 
it is true for all n < N. The function ig'^_i contains no ei(t) or e 2 (£) component. Since the 
g n are bounded, / Y,f=i 9j SW-j-i is i n L l (IR) since / is. The only way this function can 
be in L 1 is if it's exit) and e 2 (t) components are equal. This implies the result for g N , and 
the induction can proceed. 

The second statement follows from the first because the absolute values of (ei(t) + ^(t)) 
and Cj{t) for j > 3 are dominated by (1 + t 2 )^ 1 which has integral tt. 

The third statement follows since h n is obtained from g n by removal of the n th order pole 
terms. | 

We now examine G n more closely. We note that highest order pole terms in g n at t — ±i 
satisfy the recurrence relation 

Gt +1 (t) = (i^*-® ± E G f{t)G^{t) 

with Gt(t) = ± and Gt(t) = ^— r-r. From this it follows that 

1 w 4 t Ti 4 (tTi) 2 

G n {t) = i ln (e 2 „_i(t) + e 2n {t)) , (3.9) 

where 7 n satisfies the real numerical recurrence relation 

Y n—X 

7„+i = n 7 n - - E li 7n-j, (3-10) 

4 3=1 

with 7x = 7 2 = 1/4. 

By Lemma EIH the quantity (3 n = 7 n /(n — 1)! is bounded. It satisfies 

Ai+i = A " 7 E — ~~ — , — — /3n-i> (3-11) 

4 j=1 n. 
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with f3\ = fa = 1/4. From this relation and Lemma 13.21 it follows that f3 n has a limit 
(3* as n tends to infinity. To see this, suppose that the sequence (3 m is positive and strictly 
decreasing for m < n. Then, 

.<g '"-»r-" P ^<lA' 

Thus, iterating and using /3i = Pi = 1/4, we have 

2 Pf a 2(51 j2 1 

3=2 



+ 3 n(n-l) P2 3 ^ 2 j(j - 1 



- K)>a 

Therefore, the sequence {/3 n } is positive, strictly decreasing and bounded below. 
Similarly, for some constant C > and any p > 0, we have 

Pn - /W P < C — J— - V (3.14) 
Y n — 1 n + p — 1 y 

so that, 

& = /T (1 + Oin- 1 )). (3.15) 
Remark Later, we will see that P* 



iry/2' 

4 Optimal Truncation 

We begin this section by studying 

Ue,t) =ie^- H(t)^ (4.1) 

n 

where ip is given by (|2.2j) with g(e, t) = E 9j(t) e3 - We ultimately choose n = [1/e] — 1, 

j'=i 

where [A;] denotes the greatest integer less than or equal to k. 

By explicit calculation, ( n (e, t) equals e %t l^ ^'=i fj( s ) e3 ) ds $ 2 (t) times 

n— 1 2n+l n 

i^ 1 (? B + i^ 1 ii; + ie» tl /Effi^+ E ^ fc / E //,///■ > .• (4-2) 

Lemma 4.1 T/ie first term in fl^.ifl ) satisfies 

e n+1 \\G'J = 2p n e n+1 (n\). (4.3) 
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When n is chosen to be n = [1/e] — 1, the norm of the remaining terms in fl^.ifl ) satisfies 
\( n (e,t) -te n+1 G' n \\ < Ce n+1 ((n-l)\)log(n-2). (4.4) 
for some C. Thus, as e tends to zero with n = [1/e] — 1, £ n (e, t) is asymptotic to 



Proof The result ()4.3j) was proven at the end of the previous section. 
By Lemma 13.41 the second term in ()4.2|) satisfies 

e rt+1 \\h' n \\ < C'e n+1 {n-l)\ log(n - 2). 

By Lemmas 13. 1| 13 .2\ and I3.3| the third term satisfies 



,-n+l 



n-1 

/ E ft 9n-j 

i=i 



< 6^ +1 ^(n-2)! = ^e" +1 (n-2)!. 



(4.5) 



(4.6) 



We now prove that 

2n+l 



E e / E 9j9k-j-i 

k=n+2 j=k—n—l 



< C"e n+1 



(n-1)! 



(4.7) 



for some C" when n = [1/e] — 1. 

We begin the proof of ()4.7j) with a technical lemma: 



m—l 



Lemma 4.2 For positive integers I < m/ 2, V] (p!) ((m— p)!) < (m— 2/+1) (/!) ((m— 

P =i 

Proof 

m— / 

2 (pO (O-p)0 

= (/!)((m - /)!) + ((/ + l)!)((m - I - 1)!) + • • • 
+ (( m _ / _ !)!)((/ + 1)!) + (( m _ /)!)(/!). 

The first and last terms are the largest terms in this sum, and they equal (/!) ((m — /)!). The 
lemma follows since there are (m — 21 + 1) terms in the sum. | 

We apply Lemma [4.21 with p — j — 1, m — k — 3, and I = k — n — 2, along with Lemmas 

n 

13.11 13.21 and 13.31 to see that the norm of the term i e k f ^2 9j 9k-j-i in f|4.2|) is bounded 



j=k—n—l 



m—l 



by e fc - Y, P l (( m -P)0 < e fe (2n-A; + 2)((A;-n-2)!)((n-l)!)/2. 
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To prove ([4.7)1 . we now sum this quantity over k: 



2n+l 

Y e fc (2n-fc + 2)((fc-n-2)!)((n-l)!)/2 

k=n+2 

2n+1 9n h 4- 9 

= e^((n-l)!) Y 9 £ k ~ n ~ 1 ((k-n-2)\) 

k=n+2 
m=0 A 

We now fix n = [1/e] — 1 and note that this implies that 

j2 VLzJH^( m \) 



m=0 ^ 
i re— 1 

Z m=0 

1 / 12! (n- 1)! 

2 \ n n z n n ~ L 



By Stirling's formula, j\/v? < C (j/n)° e 3 \f] for some C > 0. So, the quantity on the 
right hand side here is bounded. This and (|4.8j) imply (|4.7j) . which proves the lemma. | 

We note that our choice of n = [1/e] — 1 implies that there exists a C > such that 



„ re 



< Ce. (4.9) 

i=i 

This follows from Lemmas 13.31 and 13.51 since the left hand side of ()4.9|) is bounded by 

n re re— 1 

INI ei < ttEO'-I)!^' < Tre^e^A;! < Ce. 

j=l j=l fc=0 

This allows us to define the optimal adiabatic state ^>i(e, t) associated with the eigenvalue 
-1/2 by 

Me, t) = e^/W'fe')* + ff ( c , t) S a (f) ) , (4.10) 

[l/e]-l 

where g(e, £) = 2J e3 > with the <7j(i)'s defined in Section 2. By construction, ipi(e, t) 

is normalized as t — > — oo. The optimal adiabatic state ^(e, t) associated with the eigenvalue 
1/2 is defined similarly, according to ()2.3|) . 

Me, t) = e-i^cc/W^-)* ($ 2 ( f ) + g( e , *) ) , (4.11) 

[l/e]-l 

where g(e, t) — Y, 9j(f) ^ ■ Since the entire analysis of the g/s above does not depend 
on the sign of t, it holds for the (jj(t) as well. See Remark 5 of Section 2. 
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5 Analysis of Some Integrals 



The main goal of this section is to analyze the two (quite different) integrals 
ds, 



'l±is) 



where e > is a small parameter, and m = [1/e]. By taking conjugates, we obtain the 
analogous results for 



e — is/e 

[IT is\ 



ds. 



Remark As we see below, the reason the two integrals have different behavior is that in 
one case, there is a cancellation of rapidly oscillating phases, while in the other, the phases 
reinforce one another. 

Lemma 5.1 For small e > 0, m = [1/e] and any 7 e (1/2, 1), we have 



and 

r t e is / e 



1-is) 



ds = 0(m~ 7 ). (5.2) 



Proof We begin with some preliminary estimates that apply to both integrals. We note 
that for real s, 



e is/e 



(IT is)' 



[l + s 



2\-m/2 



/oo 
ds/s m = l/(m — 1), we make an 0(m~ l ) 

error in each of the integrals, if we drop the contributions from \s\ > 1. 

Next, let 5 > be small and a = m 5 j y/m. If a < \s\ < 1, then (l+s 2 )~ m/2 < (l+a 2 )~ m / 2 . 
Thus 

f 1 (1 + s 2 )- m/2 ds < (1 - a) (1 + a 2 )- m/2 < (l + a 2 )" m/2 . (5.3) 

J a 

Now, if S is small enough, ()5.3|) is of the order of 



e -% H iW) = e -^+o(i/^-«) = (e~^) « 0(m^). (5.4) 
Therefore, if we let A m = {s : \s\ < a = m 5 /y/m}, then 

rfs = / XA m (s) + ^(m- 1 ). (5.5) 



1 ±is) m J-a (l±is) 
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We now concentrate on (j5.1|) . 

I -is \ m/2 



:i+ts)- m = n + s 2 )- m/2 



1 + is 



We separately examine the logs of the factors on the right hand side of this equation. 
First, 

log((l + *T m/2 ) = -y log(l + S 2 ). 

On the support of XA m , this equals 



™ (V - s 4 /2 + ...) = -ms 2 /2 + O^ 1 - 45 )). 



2 

Exponentiating, we have 

(l + s 2 )- m/2 = e" ms2/2 + 0(m-^-^) on the support of X A m - (5.6) 
Second, 

( / \ — is \ m / 2 \ 
log I — ) = - im arctan(s) = — i m (s + 0(s 3 )). 



l + is 

On the support of XA m , this equals 

- i m s + 0(m- (1/2 - 35) ). 

Exponentiating, we have 

/l - is \ m/2 
\l+is 

For m = [1/e] , we have 

e is/e = e ims e is(l/e-m) = jms ^ + 0{m~^ 2 ^)) = t lmS + O (m~ {1/2 ~ 5) ) , (5.8) 



v m/2 

J = e" ims + 0(m- {1/2 - 35) ) on the support of X A m - (5.7) 



on the support of XA m - 

The support of the integrand in ()5.5|) has length (^(m^ 1 / 2-5 )). Using this, ()5.6|) . (j5.7jl . 
and ()5.8|) in ()5.5|) we see that for S small enough, 

ds = / ^ , ■ ^ XA m (g)^g + 0{m- v ) 



1 + is) rn J-a (1 +is 

ft 



fl/2 „ e— 2 / 2 + 0(m-^)). 



By simple estimates on the tail of the Gaussian, this equals 

f e- ms2 ' 2 ds + 0(m-( 1 - 4 ' 5 )), 
which implies (15. 
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We now turn to (|5.2|) . The analysis is very similar, except that we use the conjugate of 
(I5.7J1 . This leads us to 



(1-is) 



ds 



-1\ 



a (1 — is) 



m-U/2-5) 



XA m {s)ds + 0{m 



ms2/2 e 2ims X A m (s)ds + Oim^ 1 -^). (5.9) 



Integrating by parts, we have 

/ e - ms2/2 e 2ims ds = — ^— e~ ms2/2 e 2ims + — [ s e~ ms2/2 e 2ims ds. 
J 2im 2i J 

Since / 

J — oc 

£, and (|5.2|) follows for 5 small enough. 



2 2 

I si e~ ms I 2 ds = — , the integral on the right hand side of (15.91) is 0(1 /m) for all 

oo m 



In the following section, we will also need the following trivial estimate for any m > 2: 
ds < C. (5.10) 





1 


/ — oo 


(1 ± is) m 



6 Proof of the Theorem 



The perturbation expansions of Section |21 generate formal approximate solutions to (jl.lj) 
associated with the energy levels =Fl/2. When truncated at n — [l/e] — 1, the expansions 
define the optimal adiabatic states ipi(e, t) and ^(e, t) by (|4.1U|) and (|4.11j) . 

As we see below, each of these optimal adiabatic states agrees with an exact solution up 
to an exponentially small error. Since ipi(e, t) and ip2{t, t) are asymptotically normalized 
for t — > — oo, it follows that they are normalized up to the same exponentially small error 
for all t. Similarly, the inner product (ipi(e, t), ^(e, t) ) is exponentially small for all t. 

We compute the asymptotic leading term to the solution of (jl.lj) that coincides with 
ipi (e, t) as t — » — oo, using the optimal adiabatic states ipi(e, t) and ^(e, t) as a basis. This 
approximation is accurate up to an error of order e~ 1//e e M , for < /i < 1/2, uniformly for 
t E R. 



We now define £(e, t) by (j4.1j) with n = [1/e] — 1. We denote the unitary propagator 
associated with the Schrodinger equation (jl.lj) by U e (ti, t 2 ). There is an exact solution 
\&i(e, t) to p.lj) that is asymptotic to ipi(e, t) as t — > — oo. We compute 

*i(e, t) - V>i(e, t) = lim U e (t, r)^i(e, r) - ^(e, t). 

r— >— oo 

In this expression, we replace U e (t, r) tpi(e, r) with the integral of its derivative with respect 
to r (with the proper constant of integration) to obtain 

*i(e, t) - <Me, t) = Jm^ - jf* £/ e (t, s) ^ if(s) Vi(e, s) + |^^i(e, s) J 
= " T ^e(*, s) C(e, s)ds. 

6 J — oo 
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To bound the difference between \l/i(e, t) and ipi(e, t), we use the unitarity of the evolution 
operator together with the expression we derived for £(e, s) in Section 0J 
This yields 

*!(e, t) - Vi(e, t) 

= - e n f e ir/{2e) J-~ {f{s) ^ 9i{s)ei)ds (G' n (r) + F n {r)) U e (t, r)$ 2 (r)rfr 



oo 



t 

e n / e" 



- e n f e ir/(2e) e £» (/(s) E "=> F n (r) tf e (t, r) $ 2 (r) dr, (6.2) 

J — oo 

where -F n (r) is — i/e n+1 times the expression (|4.2|) with the first term removed. 

We now examine (|6.1|) . which, as we shall see, dominates (|6.2|) . Using (|5.1Uj) . and (|4.9j) 

we first see that 

e / r (/W YT j= i9j{s)^)ds = 1 + 0( - e ^ 

where the 0(e) is uniform in r. 
We also have 

U e (t, r) $ 2 (r) = e -*(*-«-)/( 2e ) $ 2 (t) + 0(e) = e fr ^ ^ 2 (e, t) + 0(e), 

where the 0(e) is uniform in t and r. Using these estimates together with Proposition 13.11 
and (|5.10|) . we see that (j6.1|) equals 

- e n ^ 2 (e, t) f e ir/e G' n {r) dr + 0(e n+1 n!), (6.3) 

J — oo 

with the remainder estimate uniform in t. By (13. 9J) . (|3.10j) . ()3.11|) . ()3.15|h our choice of 
n~l/e, and Lemma f5. II with m = n + 1, we see that for 1/2 < 7 < 1, ()6.3|) equals 



+ 0(e n+l n\). (6.4) 
We deal with the prefactor before the curly bracket using Stirling's Formula 



k\ = V2Vk k k e' k (l + 0(l/fc)) 
(applied with k — n + 1) and with the easily checked asymptotics 



\LA\ = 1 + Q(l/x) as x -> 00. 



By the mean value theorem, the erf function satisfies 

erf (v^et) = erf (J[xj r) + ^ e~ r2y [x — [x]), where y G ( [x], x), r E M, 

Zy 

= erf (\/[:£]t) + 0(l/x), uniformly in r. (6-5) 
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Finally, the remainder term in ()6.3|) is of order 0(y/ee x / e ) so that the quantity (j6.1|) equals 

- P*ne-^ jerf^^tj + 1 + 0(e") J ^(e,t) + 0{e-^y/e), (6.6) 

uniformly in t for any < fi < 1/2. 

Since F n £ L , a similar analysis of ()6.2|) shows that it is uniformly bounded in t and that 
its absolute value is smaller than the estimate of the absolute value of the integral in ()6.3|) 
by a fa<;tor f C ^ M» ' 2) ^ ^ ^ ^ ^ ^ logfr - 2) ^ 

n n 

Lemma f4. II and a factor of \Jn + 1 that appears in the denominator of (|5.1jl . Therefore, by 

estimates similar to those above, ()6.2|) is of order 0(e~ 1//e ^/e ln(l/e)), as e — > 0. Hence, for 
any < /i < 1/2, 

= ViM) - Fire' 1 /* jerf^^tj + 1 + 0(e") J ^ 2 (e,t) 

+ OCe-^v^Ml/e)), (6.7) 
uniformly for t & M. 

A similar analysis of the solution ^C 6 ) of (|1-1|) that coincides with ^(cj t) as i — > — oo 
leads to the estimate 

* 2 (e,*) = ^M) + /3*7re- x / e j er f^^tj + 1 + 0(e") J Vi(M) 

+ ©(e-^v^MVe)), (6.8) 
uniformly for t £ M. 

For any < fi < 1/2, the following estimates follow directly from ([6.7)1 . ()6.8|) and the 
unitarity of U t (t,r): 

^•(e, t) = ^-(e, t) + 0( e -^), 

||Vj(e, *)|| = 1 + Oie-V't?), and 

(Vi(e,t), ^(e,t)> = 0(e- 2 M 
for j = 1, 2. Thus, by (ftJOjl and flUD}, we have 



Km ^(e,t)>| 



e 



(-l)i+i/^/(«) fl (e,(-iy+i a )<to 



1 + 0(e- 1 / £ e /1 ). 



We prove that the value of /3* is l/(7rv2) by comparing the asymptotics above as t — ► oo 
with the transition amplitude (jl.9|) . 
This proves Theorem 11.11 with 

Xj(*,t)=ip j (e,t/6) (6.9) 

when we return to the original notation e ^ e/(E S) and original time variable t*—*t/5. H 
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